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Ising model on the Apollonian network with node dependent interactions 
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This work considers an Ising model on the ApoUonian network, where the exchange constant 
~ ^/i^i^j)'^ between two neighboring spins is a function of the degree k of both spins. 

Using the exact geometrical construction rule for the network, the thermodynamical and magnetic 
properties are evaluated by iterating a system of discrete maps that allows for very precise results in 
the thermodynamic limit. The results can be compared to the predictions of a general framework for 
spins models on scale-free networks, where the node distribution P{k) ~ k with node dependent 
interacting constants. We observe that, by increasing fi, the critical behavior of the model changes, 
from a phase transition at T = oo for a uniform system = 0), to a T = phase transition 
when ^ — 1: in the thermodynamic limit, the system shows no exactly critical behavior at a 
finite temperature. The magnetization and magnetic susceptibility are found to present non-critical 
scaling properties. 
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with a different value 7 . This makes it possible, for 
instance, to devise models in the standard BA network 
that undergo finite temperature phase transition. An 
analytic expression 



I. INTRODUCTION 

Magnetic models on complex networks have quite dis- 
tinct behavior from those on regular lattices Their 
properties are of far greater importance than just a math- 
ematical curiosity, since they establish landmarks for the 
behavior of many systems, like social, economic, and 
communication networks. For such systems, the un- 
derstanding of the conditions leading to a phase tran- 
sition, or a sudden collective change in the behavior of 
the agents, is of utmost importance to avoid a breakdown 
of social structures or collective current day technological 
facilities 0, H|. 

The absence of a finite temperature phase transition 
in the thermodynamic limit T — > 00, where N is the 
number of nodes 0, Q , concomitant with the presence 
of a finite degree of magnetic ordering, stays among the 
first results that have been obtained for Ising models on 
the standard Barabasi- Albert (BA) scale- free network Q , 
where the exponent of the node distribution P{k) ~ k~'^ 
assumes the value 7 = 3. It was also observed that finite 
temperature critical behavior is found when 7 £ (3,5], 
while, for 7 > 5 the critical behavior collapses at T = 0. 
Later, an interesting interplay between critical behavior 
and node dependent interaction constants has been ev- 
idenced 0, Q- if the strength of interactions in a BA 
network, with a given value 7, is non-uniformly reduced 
according to 



7' = (7 - - m) 



(2) 



jo/{hk,r, 



(1) 



where kg is the degree of node £, the critical behavior 
moves into the universality class of the uniform model 



has been derived based on scaling arguments but, al- 
though it has been numerically verified for BA networks, 
it is not known whether its validity extends to other net- 
works. 

The purpose of this work is to investigate the effect of 
a node dependent coupling constant on the properties of 
an Ising model on the Apollonian network (AN) ■ 
This network has very special features (Til . [T3 | , includ- 
ing presenting a power law distribution of node degrees, 
with exponent 7 ~ 2.58. Previous results for Ising mod- 
els on the AN have shown that, for a variety of situations 
where both ferro- and antiferromagnetic interactions are 
allowed, phase transition in the thermodynamic limit oc- 
cur only at T = cx) AN's are constructed accord- 
ing to precise geometrical rules, which lead to exact self 
similar patterns and scaling properties. They are also 
amenable to mathematical analysis based on renormal- 
ization or inflation methods, as the transfer matrix (TM) 
formalism we will use here, which allow for the evaluation 
of its properties in the thermodynamic limit. 

These facts turn this model particularly suited for test- 
ing the existence of a finite temperature phase by mod- 
ulating the coupling constants according to Eq. ([T]). On 
the other hand, since the AN geometric rules lead to a 
well defined value of 7, there is no general free parameter 
we can use in the study to verify the validity of Eq. 
Further, it must be stressed that, despite the fact that 
AN displays power law distribution of node degree, it dif- 
fers substantially from BA network with respect to other 
topological properties, as the existence of many closed 
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loops. This is expressed, among other measures, by the 
clustering coefficient C, which is very high 0.85) for 
AN and very small (~ l/N) for the BA 

The rest of this paper is organized as follows: Section 
2 introduces the basic properties of AN networks and of 
the proposed model; details of the used TM scheme to 
evaluate the thermodynamical properties are discussed 
in Section 3. We discuss our main results in Section 4, 
emphasizing the emergence of a cross-over in /i. Finally, 
Section 5 closes the paper with our concluding remarks. 



II. APOLLONIAN NETWORK AND MODEL 

AN's have been recently introduced in the complex net- 
work framework although the original concepts 
can be traced back to ancient Greece, where the prob- 
lem of optimally filling two and three dimensional spaces 
with circles and spheres has been studied by ApoUonius of 
Perga • The complex solution to this problem, which 
amounts to placing tangent structures with well defined 
radii at precise centers, suggests the far simpler problem 
of constructing the AN. In this case, one just has to put 
a node in each circle center, and a network edge between 
the centers of each pair of tangent circles. This process 
can be followed in a recursive way in terms of the gener- 
ation g in which new circles are added to the structure. 
In this work, we consider that, at the zeroth generation 
5 = 0, three tangent circles with the same radius occupy 
the centers of an equilateral triangle (see Fig. 1). For the 
g + 1-th generation, the network construction consists in 
putting a node within each triangle of the g-th genera- 
tion, and connecting it to each of the triangle corners. 
It is a simple matter to verify that the number of net- 
work nodes N{g) and edges B{g) increase according to, 
respectively, N{g) = {33 + 5)/2 and B{g) = (39+i +3)/2. 
The average number of neighbors per node equals 6, since 
B{g)/N{g) ^ 3 in the limit g ^ oo. 

For a given generation g > 0, the largest node degree 
is kc{g) = 3 X 2^^^, where the subscript c indicates that 
such node occupies the central network position. The 
second largest degree nodes, with ke{g) = 2^-1-1, occupy 
the external corners. At any generation g, there will be 
nodes with degree k = kc{'g),'g = 1, . . . ,g and fee- The 
degree dependent node multiplicity is m{'g) = 3^^^ for 
the internal nodes, and m{g) ~ 3 Vg, for the nodes at 
outer network corners. 

As already quoted, the resulting AN is scale free. How- 
ever, it also has other properties that are typical for other 
complex network classes, as being small world (mean 
minimal path {£) ^ hiiV), hierarchical (the clustering 
coefficient of individual nodes c{k) has a power law de- 
pendence on k), and having a large clustering coefficient 
C. Because of this, systematic network clustering analy- 
sis based on several independent measures [l^ shows that 
AN does not belong to the same class as the most stud- 
ied network sets, generated by the algorithms proposed 
by Watts and Strogatz [l^ and Barabasi and Albert f6l|. 




FIG. 1: Geometrical construction of the first three generations 
{g = 0, 1, and 2) of the AN. Nodes are numbered according 
to the scheme used in [l^ 

We consider the Ising model with spins <Ti ~ ±1 placed 
on each site of the Apollonian network. Pairs of spins 
which are neighbors on the network, interact with 
coupling constants Jij. Thus, the Hamiltonian for the 
system can be written as 

Hg = Ji^j a,aj -hJ2^^: (3) 

where Jij is given by Eq.(IT]). In our previous studies, 
we have considered inhomogcneous models, in which the 
constants Jij depend on the generation g at which the 
edge, hence the second spin in the pair, was introduced 
into the network. Due to the fact that, at each genera- 
tion, the newly introduced nodes are connected to nodes 
that were introduced in previous generations, the scheme 
introduced in Ref. [l^ does not assign the values of Jij 
according to the rule of Eq. ((T]). In the following Section 
we discuss how to implement the interaction constants 
of Eq. Il]) in connection with the TM method used to 
evaluate the model properties. 

III. TM RECURRENCE MAPS 

The basic steps to implement the TM method we use 
to evaluate the thermodynamic properties have been pre- 
sented, with some detail, in one of our previous works 
p^ . However, the method needs to be adapted to the 
specific situation introduced by the more complex inter- 
action given in Eq. ([1]). Thus, let us briefly recall that the 
TM scheme amounts to write down the partition func- 
tion Z{T, h, N{g)), where H denotes the magnetic field, 
for any value of g in terms of a TM that describes the in- 
teractions between any two of the outer AN sites. In this 
process, it is necessary to perform a partial trace over all 
interaction dependent configurations. Due to the exact 
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geometric AN construction rule, it is possible to express 
the TM matrix elements at generation g + I in terms 
of the corresponding elements at generation g. In this 
framework, we basically work with a set of 2 x 2 square 
matrices 



c d 



and a set of non-square auxiliary matrices 



Pg Ig ^g ^g 
t„ u„ v„ w„ 



(4) 



(5) 



which explicitly include the dependence of the third outer 
node spin variable. As the Lg matrix elements arc num- 
bered according to the lexicographic order, the following 



relations hold: a. 



g=Pg + %: 



^g ~^ ^.9' 



dg = Vg + Wg. For more symmetrical models, and field 
independent situations, the number of independent vari- 
ables can be reduced. 

For the homogeneous systems, it was possible to write 
down a single set of recurrence relations between matrix 
elements in successive generations. Although the basic 
idea of the method remains the same, for the current 
model, it is necessary to track the way the nodes are 
reconnected when they go from 17 to g -1-1. This influences 
the change in their degrees, so that the same node will 
contribute differently for distinct values of g. 

We start the discussion of the changes in the TM 
scheme by pointing out that, besides knowing the set 
of node degrees and corresponding degree (node) multi- 
plicity, it is necessary to go one step further, and identify 
each of the P{g) = {g^ — g + 2)/2 different triangles in 
which the g network can be disassembled. In this respect, 
each triangle is characterized by the node degrees ki , kj , 
and ki of the nodes and i, respectively. For any tri- 
angle and any g > 1, there is always (only) one node with 
k = 3. Note that P{g) grows only with the square of g, 
so that, even for a complex interaction structure, there 
is practically no constraint to numerically compute these 
matrix elements for very large values of g. 

Once this set has been identified, we evaluate the 
model properties at generation g by computing the con- 
tribution to the partition function from each of these 
P{g) triangles, storing them in corresponding TM's Lg, 
i^l,...,Pig). 

To proceed further, we must consider that the g eval- 
uation is equivalent to the g — 1 one, provided we start 
with triangular units with a fourth node added at the 
central position. This way, it is possible to compute the 
contribution of the new P{g— 1) triangles, by performing 
partial trace over the contributions from the central node 
of each of these structures. The new form of the general 
recurrence relations for the matrix elements. 



{Lg+i)Zk - y^XLg)lii{L 



g)i,tk{L\)k,jti 



(6) 




is quite similar to that of the uniform model. The differ- 
ence refers to the superscripts 77, e, and k, which identify 



FIG. 2: Dependence of Jo on fi, as a result from fixing the 
ground state energy per spin uo ~ —3. 



which three TM's (corresponding triangles) have been 
put together. The same arguments can be used again, 
until we obtain one single TM that accounts for the con- 
tributions of all network nodes. 

The results we present in the next Section consider 
g < 50 which, for the largest value, is roughly of the order 
of magnitude of the Avogrado number. The adaptation 
of the uniform TM procedure to take into account the 
node dependent interaction constant depends basically 
in the identification of the basic triangular units and the 
assembling rules that combine them when one goes from 
g to g + 1. A summary of the implementation of the 
details is provided in the Appendix. 

Finally, it is important to note that the map iterations 
can be more conveniently performed if we rewrite the 
set of recurrence maps given by Eq. ^ in terms of the 
free energy and the ratio of the Lg matrix elements to 
the largest one {qg). Indeed, this avoids numerical diver- 
gences, in the low temperature region, when g increases, 
as conveniently discussed in Ref . ^ISi] . 



IV. RESULTS 

According to the previous Section, we present re- 
sults for fixed number g of generations, usually g = 
10,20,30,40, and 50, corresponding to networks with 
2.953 X 10^,1.743 x 10^ 1.030 x 101^,6.079 x lO^^ and 
3.590 X 10^'^ sites, respectively. The precise numeri- 
cal evaluation of the free energy f(T,h) allows to ob- 
tain the entropy s{T), specific heat c(T), magnetization 
m{T,h = 0), and susceptibihty x(T, = 0). It is also 
possible to calculate the ratio A1/A2 of the two TM eigen- 
values. For models on Euclidian lattice, as well as on hier- 
archical and several fractal structures [l3], this quantity 
is directly related to the correlation length ^. In the case 
of complex networks and, in particular, of the AN, the 
connection between these quantities is not so obvious, as 
the distance between two outer nodes remains always 1 



4 



C 0.4 




8 10 12 



80 
60 
40 
20 

10' 
,10° 

h^io-' 

^^10"' 

10-^ 



(a) 





(b) 



10 20 30 40 



50 



FIG. 3: Plots of the specific heat c for g = 40 (symbols) and 
g = 50 (solid line) for different values of /x: (squares), 0.5 
(diamonds), 1.0 (up-triangles) and 3.0 (down-triangles). The 
superposition of curves and symbols for g = 40 and g = 50 
indicates that c converged to its value in the thermodynamic 
limit. 



FIG. 4: Dependence of the temperatures Td (symbols) and 
Tm (solid lines) on the network generation g. (a) /i — 0.0 
(squares), 0.5 (up-triangles) and 0.8 (down-triangies). (b) 
fj, = 1.0 (circles) and 1.2 (diamonds). 



eration (7, and further on whether 0</.t<l,/i=l and 
^ > 1. So it is adequate to discuss them separately. 



in any generation. Therefore, we will discuss the behav- 
ior of l/ln{Xi/ X2), although we refrain ourselves from 
calling it ^. 

Eq. ([T|) indicates that the coupling constants Jij 
linearly depend on Jq. According to Section 2, the 
number of connections in the AN at generation g is 
Lg = {33+^ + 3)/2. If we take Jq = 1 when ^ 0, 
the free energy per spin in the g 00 limit is uq =^ —3 
at T = 0. If we fix Jo = 1 and let fi increase, the value 
of Uq decreases and, besides that, all thermodynamic ef- 
fects will occur at a lower value of T. Thus, to avoid 
choosing an adequate temperature scale to work with at 
each value of /z, we find it more convenient to choose a 
fi dependent value Jo{n), by requiring that uq = — 3,V/i. 
In Fig. 2, we show the dependence of Jo(m) on /x, which 
shows that Jq ^ exp(/i). As a consequence of this choice, 
all maxima of the specific heat occur roughly at the same 
value of T. 

Fig. 3 shows, for g — iO and 50 and /i = 0, 0.5, 1, and 
3, that the curves for c are rather insensitive to the values 
of g and /x. Moreover, they are completely smooth, with 
a Schottky like maximum at a temperature Tg . This con- 
stitutes a main difference to the results for the 7 = 3 BA 
networks [Ti]. There is reported the presence of a finite 
critical temperature, identified by a jump in the specific 
heat when > 0.5, which changes into having a diverg- 
ing slope when 0.5 > fi > 0.33. According to Eq. ^ and 
to the AN known value of 7 ~ 2.58, similar critical be- 
haviors should emerge for ^ > 0.61 and 0.61 > fi > 0.47, 
if the AN were to fall within the BA universality class. 
This clearly shows that the validity of Eq. ^ can not, 
in general, be extended from the BA to other network 
classes, even if they are scale free as the AN. 

The same calculations reveal that the resulting pat- 
terns for m,x E^nd l/^n(Ai/A2) depend, first, on the gen- 



A. < < 1 

Within this parameter interval, we observe that, like 
for the uniform model, l//n(Ai/A2) numerically diverges 
(> 10^°") for a non-zero temperature Td{g), which in- 
creases linearly with g. Since g ^ IniV, Td depends 
in a logarithmic way on the system size. If we write 
Td{g) ~ A{fi)\nN, we find that A{^) decreases with /i 
(see Fig. 4). 

In Fig. 5, the behavior of the zero field magnetization 
m{T,h — 0,N{g)), which is exactly one when T — 0, 
slowly decreases when T increases. Its behavior when 
fi = 0.5, at larger values of T, is a bit more complex than 
that for /i = (Fig. 5a). There it is clear that m suffers 
a first cross-over to an exponential decay at T^, which is 
followed by a transition to a second exponential decay, 
mediated by a larger constant, at Tm{g)- The magneti- 
zation curves for different g collapse during the first and 
second regimes. The third regime will later on be in- 
terrupted again by a smoother decay. As observed with 
Td{g), Tm{g) ~ B{ii)g, with B{n) ~ A(fj,). However, Td 
and Tfn do not coincide. The second part of the magneti- 
zation curves, where m(T^O, N{g)) overlaps for different 
values of g, extends over wider T intervals when g in- 
creases. This shows that, in the thermodynamical limit 
g —>■ 00, the value of m will follow the second exponential 
decay when T 00. Nevertheless, as observed for Td, 
this region grows logarithmically with the network size. 

The behavior of x is strongly correlated with that of 
TO. It vanishes when T ^ 0, then it grows with T, shows 
a first maximum at a 17 independent Ts, and a second 
g dependent maximum at Tm- As for m, the curves for 
larger values of g overlap for much larger distances. The 
maxima of the x curves are described by an universal 
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function, as can be observed in the very precise re-scaled 
curves in Fig. 6a, which shows that the scahng exponents 
increase with ii. Note that only the value of x needs to be 
scaled by the corresponding maxima, while the location 
at the temperature axis is corrected by shifting the scale 
by Tm{g). This excludes any possibility of having a crit- 
ical phenomenon associated with susceptibility maxima. 



This value of determines a crossover in the behavior 
of the system, which is reflected both in m and x- This 
change can be noticed in Fig. 3, which shows that — 
1) = = 1) = 0, i.e., the temperatures associated with 
the maxima of the susceptibility and the divergence of 
l/ln{\i/X2) become independent of the system size. The 
precise value of Td depends, of course, on the threshold 
value of the numerical divergence. However, by plotting 
the value of l/Zn(Ai/A2 as function of l/T, we notice a 
linear dependence in the T ^ limit, suggesting that 
Td = 0. 

However, this new behavior cannot be associated with 
the emergence of criticality. First we recall that Fig. 
3 does not indicate any change in the Schottky profile 
and, second, we see that Tm > Tg. Finally, the x curves 
in the region around Tm, which shows a perfect scaling 
with respect to g with scaling exponent 1, are completely 
smooth (see Fig. 6b). Note that the horizontal axis in- 
dicates that it is not necessary to shift temperature as 
in Fig. 6a. Note that the two maxima, which were ob- 
served when /i < 1, have merged together, and that the 
large temperature side is characterized by an exponential 
decay. 

As for the previous n interval, the behavior of m is 
strongly correlated with that for x- It is characterized 
by a single exponential decay after Tm, with a very large 
constant, as shown in Fig. 5b. 



C. /i > 1 

In the last range of parameter values, Td and T,„ de- 
crease with respect to g. As shown by Fig. 4c, both val- 
ues converge exponentially to with respect to g. Fig. 
4c also shows that the rate in the exponential increases 
with /i. 

Therefore, the behavior of l/^ri(Ai/A2) is different from 
that one observed for = 1, the divergence when T — > 
becoming slower at increasing values of g. This suggests 
that, when (7 — > cx), any collective spin ordering is weaker 
than that of an Ising chain, rather typical for a param- 
agnetic situation. 

The shape of the m curves becomes completely differ- 
ent. The stable plateau at m = 1 for a finite tempera- 
ture interval, which survived until /i = 1 , disappears as g 
increases, indicating that no spontaneous magnetization 



m 



10" 
10' 

10-^ 
10-^ 

10-^ 
10-^ 
10"' 
10"' 



i 




(a) 








- 
















\ ^ ■ - 












iii \ \ 






\ \ 












iv 







10 



20 30 
T 



40 





FIG. 5: Behavior of m(T, h = 0; N{g)) against T for different 
values of g when p = 0.5, 1.0 and 1.5. (a) /i — 0.5, g — 20, 30 
and 40 indicated by dots, dashes and solid lines. Four regions 
characterized by different behavior are obtained: Two of them 
are g independent, the second of which has an exponential 
decay. The third region starts with a cross-over to a second 
exponential regime, (b) fi = 1.0, same symbols as in (a). The 
first exponential region has disappeared, (c) jj, = 1.5, g = 
10, 15 and 20 indicated by dots, dashes and solid lines. The 
first region in (a) has disappeared. Two T intervals separated 
by a g dependent crossover temperature are observed. As g 
increases, m vanishes for any T > 0. 
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FIG. 6: Main panels show scaling properties of xiT, h — 
0;N{g)) against T for different values of (? = 40 (solid), 
45(dashes), and 50(dots), when fi = 0.5 (a), 1.0 (b), and 
1.5 (c). The inserts show the same curves in the original vari- 
ables, (a) Scaling properties of x as function of T — Tm, with 
exponent 'y/u = 0.607, are observed in the region iii of the 
magnetization m (see Fig. 5a). (b) Scaling properties of x as 
function of T, with exponent 7/1^ = 1.0, valid in the region 
iii, with exponential decrease of m, as discussed in Fig. 5b. 
(c) Scaling properties of x as function of T, with exponent 
j/i/ = 1.0. Scaling is valid in the two regions shown in Fig. 
5c. The horizontal axis is in logarithmic scale to evidence 
Curie's law with g dependent constant Cg. 



exists for a r > 0. This suggests that, in the thermody- 
namic Hmit, TO = OVT > 0. Our results do not allow to 
assert that to = also for T = 0. 

The behavior of x supports the conclusion of a low tem- 
perature paramagnetic phase. It increases very rapidly 
from to a maximum value at Tm, followed by a l/T law. 
Since Tm goes exponentially to zero, a Curie law prevails 
for large g. This is shown by the scaled curves in Fig. 6c, 
which indicate that the Curie constant C depends on g. 



V. DISCUSSION AND CONCLUSIONS 

The results we obtained for the magnetic behavior of 
the Ising model with node dependent interaction con- 
stants reveal a quite rich picture, although no critical be- 
havior at a finite temperature has been identified. The 
properties of specific heat show that the g dependent 
curves converge very rapidly to a well defined value in the 
thermodynamic limit. On the other hand, magnetization 
and susceptibility indicate a much more complex behav- 
ior which, for certain temperature intervals, are heavily 
dependent on the value of g. 

The TM method allows for the comparison of to and 
X for different values of g, which leads to the identifica- 
tion that part of the results are due to finite size events. 
The curves showing such effects are amenable to very pre- 
cise collapsing by adequate scaling expressions, similar to 
critical points in magnetic models on Euclidian lattices. 
This includes the dependence of characteristic values of 
the temperature (Td and T„i) 

The behavior of the system in the region < 1 is close 
to that observed for magnetic system with uniform in- 
teractions on BA networks: only an ordered phase is ob- 
served at any value of T. /i = 1 characterizes a crossover 
in the behavior of the system, as for /i < 1 the mag- 
netization vanishes, for any value of T, when g — > 00. 
This region reveals a typical behavior of a genuine para- 
magnetic system. This pictures is corroborated by the 
behavior of Xi as one finds that a Curie law is valid in a 
limited region close to T = 0. For larger values of T, the 
decay is characterized by an exponential decay. 

In the context of complex networks our most important 
finding is that the relation of Eq.(l) between effective 
topology and interaction strength proposed in Refs. [7,8] 
does not have general validity for all scale-free networks 
since the Apollonian case behaves differently. 
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VII. APPENDIX 

As discussed in the Section II, for any generation g, 
the node occupying the central position of the AN has 
the largest degree kdg), the value of which results form 
the difference equation relating the values of kc at two 
successive generations: kc{g+ 1) = 2kc{g)- The degree of 
the nodes at the external corners obey a similar equation, 
namely: kf,{g + 1) = 2ke{g) — 1. 

The AN can be disassembled in triangles, in such a way 
that each node i of degree ki belongs to ki triangles. The 
only exception refers to the nodes at the external corners, 
which have degree ke{g) but belong to ke{g) — 1 triangles. 
Each such triangle can be characterized by the degree of 
its three nodes. The number P{g) of different triangles 
at generation g can be expressed in terms of p{g) and 
T{g), respectively the number of triangles that does not 
include (includes) an external node: P{g) = p{g) +T{g). 
Since they obey the relations p{g + 1) = p{g) + T{g) — 1 
andT(g + l) =T(g) + l, we obtain p{g) = [g'^ - ig + 2)/2 
and T{g) = g, from which the expression for P{g) antici- 
pated in Section III follows. The number p{g) can be fur- 
ther decomposed in terms of /,(<?) and respectively 
the number of different triangles that includes (does not 
include) the central node at generation g. It is a sim- 
ple matter of inspection to see that L(g) — g — 2 and 
= (s' - 55 + 6)/2, Vg > 2. 

For the purpose of computing the TM's, it is neces- 
sary to identify the distinct triangles present in the AN. 
This proceeds by the collection Tg_(;(fci, fc2, ^3), where g 
indicates the generation, C is a number e [l,P{g)\, and 
ki indicate the degrees of the nodes at the vertices of the 
triangle. Tg^Q(ki, k2, k^) are recursively defined according 
to the following rules: 

1) fcg = 3, V5, VC- ^ 

2) Since r((/) = g, one single new triangle containing 
an external node is introduced into the network, Vg. We 
use C = P{g — 1) -I- 1 to characterize it, and note further 
that ki — kf,{g), k2 — kc{g). Any such triangle retains 
the ki values for all further generations, i.e.: 

Tgxiki,k2, fcs) = T<,,c'(fci, fc2, ks), VC' > C- (7) 

3) V5, there are p{g) different triangles, among which 
K(g) have been introduced in previous generations. They 
will be characterized by the same values of ki, so that 
Eq. ([7]) also holds for this subset. The remaining L{g) 
new triangles are numbered according to the rule: C = 
P{g — !)+£,£ = 2, 3, 5 — 1. For each value of £, we 
set ki = kc{g), = kc{i). 

The final step consists in establishing the rule to com- 
bine the contributions to the partition function from 
three distinct triangles at generation g to obtain the par- 
tition function at generation g + 1 according to Eq. 
If the properties of the systems are to be computed un- 
til a chosen value g, we are required to start with P{g) 
distinct triangles, precisely identified as discussed above. 
Then, as discussed in Section III, it is necessary to define 
a map that selects the proper values of rj, e, S € [1, P{g)] 



used to perform the trace over the common central node 
of the triangle a £ [1, P{g — 1)]. So let us note that 



[1, P{g)] = [1] U 1 U [^(5 - 1) + 1, P{9 - 1) + J - 1] 

(8) 

Then, the values of 77, e, 6 are given, as function of a, 
by the following expressions: 



a 



: 1 : 

1, 6 = 5 = 2; 

-PU-l) + l,je[2,g^l] : 
2,6 = P{j) + 1,5 = P{j) + 2- 



a G [P{j - 1) + 2, P{j - 1) + J - 1] , J e [2, 5 - 1] : 
Tj = P{a- P{j - 1)) + 2, e = P{j) + 2,d = a + j. 

As the AN is self similar, these maps are also valid for 
all forthcoming partial trace operations, until only one 
single triangle is left. At this step, the remaining TM 
contains the contribution form all spin configurations, 
from which the thermodynamical properties follow. 

With the help of these relations, a set of recurrence 
maps can derived from Eq. which allow for the eval- 
uation of the free energy and its derivatives: 



Jg+1 



exp(-2/3/i)) 



l + <ga;|,Xgexp(-2/3/i)); 



V 5 { V ^ 



X 



2,g+l 



X 



3,3+1 



^-4,3+1 



l + <3^i,9^l9exp(-2/3/i)); 



^2,34,3^,3 + 4',g^3,3^7,g exp(-2/?fe)) 

l + xlgxlgxlg exp(-2/3/i)) 



cL-4,3 + ^l,g^l,g^i,g exp(-2/3/l)) 



l + xl„xl„xl„eM-'^Ph))-, 



X 



5.3+1 



,g^6,g^2,g 



''5,3-^7,3-^3,3 ' 



-2f3h)) 



6,3+1 



l + <,^lX,exp(~2/3/i)); 



A,g^Xg^i,g + 3^7,3^5,3^5,9 exp(-2/?fe)) 

l + <3^i,X9exp(-2/3/i)); 



(9) 



(10) 



(11) 



(12) 



(13) 



(14) 



(15) 



8 



'7,3+1 



■ Xr^gX^gX^g cxp(-2/3/i)) 



l + XlgXlgXl^eM-m)). 



(16) 



Pg Pg Pg 



(18) 



In the above relations, the following variables have 
been used: 



X7,g = 



Wo 



1g, 



Pg 



' 9 . 



Xl^n — j X2 n — ^ X^ q — J 
Pg Pg Pg 



(17) 



(19) 
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